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Abstract 



r~| I We rigorously study the long time dynamics of solitary wave solutions 

of the nonlinear Schrodinger equation in time- dependent external poten- 
r~| . tials. To set the stage, we first establish the well-posedness of the Cauchy 

"j^ I problem for a generalized nonautonomous nonlinear Schrodinger equation. 

We then show that in the space- adiabatic regime where the external poten- 
tial varies slowly in space compared to the size of the soliton, the dynamics 
^ I of the center of the soliton is described by Hamilton's equations, plus terms 

^> ' due to radiation damping. We finally remark on two physical applications 

■ of our analysis. The first is adiabatic transportation of solitons, and the 

^ i second is Mathieu instability of trapped solitons due to time-periodic per- 

turbations. 

o . 

^ 1 Introduction 



1.1 Heuristic discussion and overview of earlier results 

In the last few years, there has been substantial progress in rigorously under- 
standing solitary wave dynamics of the nonlinear Schrodinger equation in time- 
independent external potentials. The basic picture is that in the semi-classical 
limit, the dynamics of the center of the soliton is described by Hamilton's (or 
Newton's) equations, plus terms due to radiation damping. This is a beautiful 
example where the solitary wave solution of the nonlinear Schrodinger equation 
behaves like a classical point particle in a suitable limit. Rigorous confirmation of 
this picture in time-independent potentials has been given in [U [2] for the Hartree 
equation, and in [3] for local nonlinearities. The case of general nonlinearities has 
been studied in [1] ; see also O [6l [7] . 

In this paper, we rigorously study the dynamics of solitary wave solutions 
of the nonlinear Schrodinger equation in time- dependent external potentials. As 
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far as we know, this problem has not been studied in the hterature so far. We 
show that when the potential varies slowly in space compared to the size of the 
soliton, which we call the space- adiabatic regime, the center of mass motion of the 
soliton is almost like that of a classical point particle in the external potential, 
independent of the rate of change of the potential with time. We also show that 
this picture holds for much longer time scales (0(| logh\/h), h ^ 1,) than the 
one given in [1] {0{h~^), h <^ 1). Along the way, we discuss sufficient conditions 
for the well-posedness of a generalized nonautonomous nonlinear Schrodinger 
equation with time- dependent potentials and nonlinearities. We finally sketch 
two physical applications of our analysis. The first is adiabatic transportation 
of solitons, and the second Mathieu instability of trapped solitons due to time- 
periodic perturbations. Our analysis relies on important developments in the 
nonlinear Schrodinger equation during the past two decades, [H [H [TOl [TTl [HI [151 
[ISl[IZl[ISl[ia[2ni[2Il[221[2a[2Sl[31[2], particularly [1]; see also for ^ [2H] for 
comprehensive reviews. 

We note that our analysis regarding the well-posedness of the generalized 
nonautonomous nonlinear Schrodinger equation is of some independent inter- 
est, especially that rigorous investigation of such equations is poor compared to 
the autonomous one, despite the former's relevance to very many experiments 
in quantum optics and Bose-Einstein condensates, where experimentalists can 
change various parameters with time. 

1.2 Notation 

In the following, L^{I) denotes the standard Lebesgue space, 1 < p < oo, with 
norm 

\\f\\Lv = {jdx\f{x)\^)K feLPiI),p<^, \\f\\L-=esssup{\f\), feL^il). 
We also define 

||/||lp(/,L9(j)) := II WfWLiiJ) ||lp{7)- 

For 1 < p < oo, is the conjugate of p, i.e.,l/p + 1/p' = 1. We denote by (■, ■) 
the scalar product in L^(M^), 

{u,v) = Re uv, m,wgL2(R^), 

as well as its extension by duality to Y x Y', where Y and Y' are complete metric 
spaces such that Y ^ ^ Y', with dense embedding 

^Unlike the case when the external potential is time- independent, the semi-classical limit 
and the space-adiabatic limit are not equivalent when the potential is time-dependent. Instead, 
the semi-classical limit, which is equivalent to scaling both time and space t — > ht,x hx, 
/i <C 1, is a special case of the space-adiabatic limit (where only space scales as a; — *■ hx, h <C 1). 
These statements will be made precise in the following sections. 
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Given the multi-index a = (cti, ■ ■ ■ , on) G N^, we denote |a| = Fur- 
thermore, := (9"i ■ ■ ■ (9^^. 

For 1 < p < oo and s G N, the (complex) Sobolev space is given by 

W''P{R^) := {u G 5'(R^) : G ^^(M'^), |a| < s}, 

where is the space of tempered distributions. We equip W^'^ with the 

norm 

\\u\\w-,P = \\d^u\\LP, 

a,\a\<s 

which makes it a Banach space. Moreover, W~^'^' is the dual of W^'^. 
We denote by 

H''P(R^) := {u G S'{R^) : J^-\l + \k\^)ij^u G ^^(M^)}, s G M, 1 < p < oo, 

where JF stands for the Fourier transform. The space H'^'P is equipped with the 
norm 

\\u\\hs,p = \\J^-\1 + \k\''f^Tu\\LP, u G //^'P(M^), 

which makes it a Banach space. We use the shorthand H'^'^ = H^. 
Given / and g real functions on M^, we denote their convolution by -k, 

f-kg{x) := / dy f{y-x)g{y). 



Given x G M^, we denote ||x|| := yXlili^i- 
1.3 Description of the problem 

In this paper, we study the long time dynamics of solitary wave solutions of 
the nonlinear Schrodinger equation in time- dependent external potentials. The 
nonlinear Schrodinger equation is of the form 

zdtij = {-A + VH{t,x))ij-f{i:), (1) 

where t/^ : M x — C, x G M.^ denotes a point in the configuration space, t G M 
is time, dt = A = J2f=i '§2~ the A^-dimensional Laplacian, Vh is the external 

potential, such that 

Vh{t,x) = V{t,hx),h G M+, 
and the nonlinearity / is a mapping on complex Sobolev spaces such that 



3 



/(O) = 0, and f{ip) = /{ip), where ■ denotes complex conjugation. Typical 
nonlinearities are local ones 

fii;) = \\i;\^^P, A > 0, < s < ^, (2) 

and Hartree (nonlocal) nonlinearites 

/W = A(iy^|^n^,A>0, 

where W is continuous, positive, spherically symmetric, and tends to zero as 
oo. Note that the above nonlinearities are self-focusing, and a general 
nonlinearity can be a sum of both local and nonlocal ones. The external potentials 
that we consider in this paper satisfy 

Vh{t, x) = V{t, hx), V{t, x) e C\R, C\R^)) (3) 

such that 

d^V e L°°(M^) + LP(M^)), |a| < 1, (4) 

p > Y,p > 1, and 

gL°°(R,L°°(M^)). (5) 

General assumptions about the model will be discussed in detail in Subsection 
12.11 We also show that the Cauchy problem with these assumptions is (globally) 
well-posed in in Section [3l 

When V = 0, the nonlinear Schrodinger equation ([T]) with nonlinearities as 
given above admits solitary wave solutions, which are stable stationary, spheri- 
cally symmetric and positive solutions 

V,{x,t) :=e^(^-(^-'^)+^)r/^(x-a), (6) 
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where a := {a, v, 7, /i}, a = vt + ao, 7 = /it-|-^t + 7o, with 70 G [0, 27r), Oq, v G 
and /i G M"*^, constant, and 77^ is a positive solution of the nonlinear eigenvalue 
problem 

(-A + ^^)v,-f{l^,) = o. (7) 

The solution stands for a solitary travelling wave with velocity v, center a and 
phase ^v-{x — a) +7, and the size of the soliton is oc in the sense that 77^ ~ 

e-v^ll4 as 11x11 ^cx), see0[IOl[IIl[Tll[T6l[I7l[l8l[T9l[2Ol[ai[23 
We consider in this paper potentials which vary slowly in space compared to the 
size of the soliton, i.e., 

Supp\W{x,t)\ 
sup — < 1, 

tm 

which corresponds to the space-adiabatic limit if we set the size of the soliton to 
0(1). 

We now state a rigorous result for the special class of nonlinearities discussed 
above. A more general result. Theorem [21 will be stated in Subsection 12.21 after 
listing general assumptions. 
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Theorem 1. Suppose the nonlinearity f is given by and that the external 
potential Vh satisfies Let Iq be any closed, bounded interval in M"*". For 

h <^1, suppose the initial condition ipo satisfies 

\\e-''2^°-^{iJo-Vao)\\m<h, 

for some (Tq G x x [0, 2tt) X Iq. Then, for small enough h <^ 1, there exists 
an absolute positive constant C, independent of h, but possibly dependent on Jq, 
such that for times < t < (7-1^^^^^, the solution to the nonlinear Schrddinger 
equation (QP with initial condition ipQ is of the form 

^(x, t) = e*(5^-(^-")+^)(r/^(x - a) + w{x - a, t)), 

where 

\\w\\hi = 0{h^), 

and where the parameters v, a, 7 and fi satisfy the differential equations 

ldtV = -{VV)it,a) + Oihl), 
dta = v + 0{hl), 
da = fi-V{t,a) + ^v^ + 0{hl), 
dtfi = 0{hl). 

In other words, for initial conditions close enough to a solitary wave solution, 
and for external potentials which vary slowly compared to the size of the soliton, 
the center of mass motion of the solitary wave is determined by Hamilton's (or 
Newton's) equations of motion for a point particle in the external potential, up 
to small corrections corresponding to radiation damping. The same result holds 
for more general nonlinearities, see Section [2l 

The organization of this paper is the following. In Section [21 we list the as- 
sumptions on the nonlinearity and the potential and state the main result of the 
paper. We also discuss models where the various assumptions are satisfied. In 
Section [3l we discuss the well-posedness of a generalized nonautonomous non- 
linear Schrodinger equation, such that ([1]) corresponds to the special case when 
only the potential is time-dependent. We then recall basic useful properties of the 
nonlinear Schrodinger equation and the soliton manifold in Section HI In Section 
we prove the main result of the paper. Theorem [21 Finally, in Section [SI we 
discuss two physical applications of our analysis. The first is the adiabatic trans- 
portation of solitons, and the second is Mathieu instability of trapped solitons 
due to time-periodic perturbations. 
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2 Main Result 

In this section, we precisely state the main result of this paper after listing our 
assumptions. 

2.1 The Model 

We now list our assumptions and discuss models where they are satisfied. 
(Al) Nonlinearity. The nonlinearity f = fi + ■ ■ ■ + fk such that 

/, G C'iH\R'', C), H-\R'', C)), J = 1, ■ ■ ■ , A;, 
satisfy the following. 

with Fj = fj, where the prime stands for the Frechet derivative. 

^'^^•^P,^^), ([2,oo],iV=l), 

such that VM > 0, 3 a finite constant Cj{M) such that 
\\f,iu)-f,{v)\\^.,<C,iM)\\u-v\\Lr,^ 

\/u,v E H^(R^X), \\u\\m + \\v\\m < M. Furthermore, 

lmfj{u)u = 

almost everywhere on ]R^,Vn G H^{R'^,C). Let F := Yl'j=i^j- -^^^ every 
M > 0, there exists a positive constant C(M) and e G (0, 1) such that 

F{u) < ^-^\\u\\hi + C(M), Vn G H\R^) 
such that II^^IIl^ < Furthermore, 

sup \\F"{u)\\B(m,H-^) < oo 

\M„i<M 

sup \\F"'{u)\\H^^B{m,H--^) < oo, 

||u||^i<M 

where B denotes the space of bounded operators. 
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(A2) External Potential. The external potential satisfies 

Vh{t, x) = V{t, hx), V{t, x) G C\R, C'^{R^)) 

such that 

d^V e L°°(M, + LP(M^)), |a| < 1, 

p > > 1, and 

dtVh{t,x)EL'^iR,L^iR'')). 

(A3) Symmetries. The nonlinearity F satisfies F{T-) = F{-), where T is a trans- 
lation 

Tf : ^u(x-a), a G R^, 

a rotation 

: ^ u{R-^x), R G 50(iV), 
a gauge transformation 

: ^ e'^M(x), 7 G [0,27r), 

a boost 

or a complex conjugation 

T"^ : u{x). 

(A4) Solitary Wave. 3/ C M such that V/i G /, the nonlinear eigenvalue problem 

(-A + /i)r/^-/(r/^) = 

has a positive, spherically symmetric solution 77^ G L^(M^) fl C^(]R^), such 
that 

|||a:|^?7^||L2 + ||knV?7^|||L2 + || |x|^9^?7^||l2 < cx),V/i G /. 
(A5) Orbital Stability. The solution r]^ appearing in assumption (A4) satisfies 

I dx r]l> 0,^11 EL 



(A6) Null Space Condition. We define 

:=-A + /i-/(r7^), 

which is the Frechet derivative of the map ip {—A + n)ip — f{ip) evaluated 
at rjfj^. For all /i G /, the null space 

A/'(£^) = span{iri^, 9^ .77^, j = 1, ■ ■ ■ , A^}. 



Remark 1. Assumptions (Al) and (A2) are sufficient to establish global well- 
posedness of the nonlinear Schrddinger equation in H^, see Sectionl^ Moreover, 
(Al) implies 

\F{u + v)-F{u) - {F'{u),v)\ < C{M)\\vfHi 

\Fiu + v)-Fiu) - {F'iu),v) - ^{F"iu)v,v)\ < C(M)||i;f^, 

\\F\u + v)- F\u) - F"{u)v\\h-^ < C(M)||^||^„ 

for any u E H'^(R'^) and v G if^(M^) such that \\u\\hi + \\v\\m < M. 
Remark 2. Assumption (Al) is satisfied, for example, if 

F{u) = \jdx G{\u\^) + ^ 

where G{r) = /J" dsg{s), such that g G C^{R+) with g{s) < C{1 + a G [0, |), 
1^3*^^(5)1 < C(l + s'?-'^), = 0, 1, 2, g G [0, iV > 3, g G [0, oo), iV = 1, 2, and 
W G LP + L'^,p > f ,p > 1, such that max{0,W) G L'{R^) + L°°{R^),r > f , > 
l,iV > 2, see E 



Assumption (A3) follows ifW{r) = W{\r\). Assumption (A4) is satisfied for 
local nonlinearities if 

— oo < \mig{s) < n 

s— ♦O 

-OO < lim s-°^(s) < C, 

s— too 

where < a < 2/{N - 2), when N > 2 and a G (0, oo) if N = 1, 2, such that 



3( > 0, such that / dsg{s) > 
Jo 



see for example fT4\ HE, HM 13 W^ - Moreover, (A4) is satisfied for nonlocal 
nonlinearities if, in addition to the above. 



N 

W G Ll^, ^ ^ Y' W as \\x\\ oo; 



see l2l\ [77], Assumption (A5) imply orbital stability of the solitary wave 

solution, see fT^ ITB^. It is satisfied for local nonlinearities f{ip) = Xlipl'^ip, s < 
Assumption (A6) is satisfied for local nonlinearities if 

g'{s)+g"{s)s'>0, 

ortfN = l, gCf. 
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2.2 Statement of the Main Result 



In this subsection, we state the main result of this paper, which will be proven 
in Section [5l 

Theorem 2. Suppose assumptions (A1)-(A7) hold. Given h> such that h <^ 
1, suppose there exists ctq = {oq, Vq, 70, /xq} ^ x x [0, 27r) x I such that the 
initial condition ipQ G with 

||e->-^(^o-r/.o)||^^i </i. 

Fix e G (0,1). Then, for h <^ 1, there exists an absolute positive constant C, 
which is independent of h and e, but which might depend on (Tq, such that the 
solution to the nonlinear Schrddinger equation (J\) with initial condition V'o can 
be written, for all time t G [0, Ce\ log h\/h), as 

ij{x, t) = e^-("-")+'^(r/^(x - a) + w{x - a, t)), 

where \\w\\h^ = 0{h^~^), and the parameters a, w,7 and fx satisfy the differential 
equations 

dta = v + 0{h^~') 

dtv = -2{VV){t, a) + 0{h'^-') 

da = fi- v{t, a) + + 0{e-') 

We will prove this theorem in Section [5l We now discuss the well-posedness 
of a generalized nonautonomous nonlinear Schrodinger equation. 

3 Well-posedness of a generalized nonautonomous 
nonlinear Schrodinger equation 

In this section, we discuss the local and global well-posedness of a generalized 
nonautonomous nonlinear Schrodinger equation with time-dependent nonlineari- 
ties and potential. We treat the potential and the nonlinearity as time-dependent 
perturbations. The application in this paper corresponds to the special case when 
only the external potential is time-dependent. 

Consider the problem corresponding to a generalized nonlinear Schrodinger 
equation 

iM = -Ai; + g{t,^), ^(t = O) = 0, (8) 

where g contains both the potential and the nonlinearity. Note that g can also 
depend on x G M^, but we drop the explicit dependence when there is no danger 
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of confusion. In what follows, we say that (g, r) is an admissible pair if 

2N 

^e[2,^^), (rG [2,oo],iV=l) 
2 1 1 

- = Nil~ -) (9) 
q 2 r 

We make the following assumptions on g. 

(Bl) The nonlinearity g = gi + ■ ■ ■ + gk such that gj G C(M.,C{H^,H~^)),j = 

(B2) There exist admissible pairs {qj,rj), j = 1, ■ ■ ■ k, such that, for every T,M > 
0, there exist a constant C(M) independent of T, and P independent of T 
and M, such that 

- < C(M)(1 + T^)||m - 

for all G i/"*^ with + ||'y||_f/i < Af, and \t\ < T, where r' is the 

conjugate of r,i. e., 1/r + 1/r' = 1. Furthermore, 

< c(M)(i + r'^)(i + ll^ll^i,.,) 

for all u G //^(M^) H Vri'^(M^) such that < M and |t| < T. 

(B3) Im5fj(t, n)n = 0, j = 1, ■ ■ ■ , A;, almost everywhere on R^, for all t G M and 
u G H\ 

(B4) There exists a functional Gj G C(M, M)) with G^- = gj, where the 

prime stands for the Frechet derivative. We let G = Gi + - ■ ■ Gk- For u G H^, 

\dMt,u)\<ci\\u\\L^m, (10) 

where G depends only on ||m||l2 and the real function / G L°°(]R) such that 
l{t) <1 for almost all t G M. 

(B5) For all M > 0, there exists G{M) > and e G (0, 1), both independent of 
t G M, such that 

\Git,u)\<^\\u\\l,+GiM), (11) 

uniformly in t G M, Vm G H^, such that ||m||l2 < M. 

We first prove local well-posedness of the Cauchy problem by extending Kato's 
method, which is based on Strichartz estimates and a fixed point argument, 
[25] . Proving global well-posedness for data which are not necessarily small 
is a little bit more delicate, since energy is not conserved. 
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Proposition 1. Suppose g satisfies assumptions (B1)-(B3). Then the following 
holds. 

(i) For every G i/^(M^), there exists a unique, strong -solution u of 
which is defined on a maximal time interval (— T*, T*), such that there exists 
a blow-up alternative, i.e., if T* < oo, ^ oo as t T*, and if 

< oo, —> oo as t \ —T^. Moreover, 

for all admissible pairs (a, b). 
(a) The charge is conserved, 

\Ht)\\L2 = 

for allte (-T,,T*). 

(Hi) u depends continuously on ■ If (pn in H^, and ifun is the maximal 
solution of l^j corresponding to the initial condition ipn, then Un u in 
C([— S*, T], L^(M^)) for every compact interval [—S,T] C (— T*,T*) and 
PG[2,^) (pG [2,oo),iV=l). 

Proof, (i). We set r = max(ri, ■ ■ ■ ,rfc) and consider the admissible pair {q,r), 
which satisfies (j9]). For fixed M,T > 0, we introduce the space 

Y := {n G L°^((-T,T),iJi(M^)) n L%i-T,T),W'''\R'')) : \\u\\Lo^a-T,nm) < M, 

\\u\\li{{-T,T),W^'^) < M} 

with distance 

d{u, v) := \\u - t'||L«=({_r,T),L2) + \\u - v\\li{{-t,t),l-)- (12) 

We note that (Y, d) is a complete metric space. To see this, consider a se- 
quence {un}neN C Y such that d{un, u) ^ as n ^ 00. Then there exist 
two subsequences {un^.} and {uny} such that Un^{t) u{t) in L^(]R^) and 
Uny{t) '"'^ u{t) in for almost all t G {-T,T). It follows that 

— ^00 

\\u\\li{{-T,T),W^''-) < linimf ll^inj^, ||L9((-r,T),VKi''-) 

and u G L°°((-T,T),i7i(M^)) n L'?((-T, T), see for example [271, 

Chapter 2. Therefore, m G Y. 

We divide the proof into five steps. 
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Step 1. Existence . Consider the selfadjoint operator Hq := —A, with domain 
i7^(M^), and let U be the unitary propagator generated by Hq, which is given by 

U{t,s) = e-'^'>^'-'\ t,seR. 

Given G H^{M.^), we introduce the mapping (f which is given by 

if{u){t) ■.= U{t,0)(f)-i [ ds U{t,s)g{s,u{s)). (13) 

Jo 

We will show that for a suitable choice of M and T, is a strict contraction on Y. 
Existence of a solution of (IHl) will then follow from Banach's fixed point theorem. 

For Tj appearing in Assumption (B2), we choose qj such that {qj,rj) are 
admissible pairs, j = 1, ■ ■ ■ , k. Applying Holder's inequality (in space) gives 



III II II .j{r-2) II II rj{r-2) 



and it follows by applying Holder's inequahty in time that 

2(r-r) r(r,-2\ 



III < W II '■j^''"^^ II II II 2) 

r'llL''J({-T,T),iy^'''i) — \\'^\\L°°{{-T,T),m)\\ W^W Li{{~T,T),W^^^)- 



Therefore, if m G Y, then 



2{r-r) r(r-2) 



\U\ 



L''m-T,T),W 



< M'-i('-2)M'-^('-2) = M, j = !,■■■ ,k, 



and u G L''^{{-T,T), W^'""^ (R^)), j = !,■■■ ,k. It follows from Assumption (B2) 
that gj G L''^{{-T,T), ly^'^^(M^)) such that 

for u G Y. It follows that 



\g,{t,u)\\^,>^^^^^^^^^,.,.^<CiM)il + T^)iT'^ +M)T^, j = l,---k. (14) 
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Together with Strichartz estimates H, we have 

||</^(M)|U^((_T,T),Hi) + ||^(w)||L.((-T,T),Tyi.O < C'U\\h^+CC{M){1 + M)T\ (15) 

for M G Y and T < 1, where C is a positive constant independent of T and M, 
and 5 = minjg{i^...fc} '^^ . Note that it follows from ([9]) that 5 > 0. Furthermore, 

it follows from Strichartz theorem that (p{u) G C{[—T,T],H^{W'^)). We choose 
M and T < 1 such that 

C"||0||j^i < M/2, C"C(M)(1 + M)T^ < M/2. (16) 

For this choice of T and M, </?(u) G Y for all u G Y. 
It also follows from Assumption (B2) that 

\\g,{t, u) - g,{t, ^)ll^,,((_^_^)^^.^.) < C(M)(1 + T^)\\u - vW^ya-W^^), 

for u,v G Y, and hence 



Applying Strichartz estimates, we have 

\\ip{u)-ip{v)\\L^i-T,T),LnH^i^)-Viy)\\L^{{-m^ 

where d{-,-) is defined in f|T2|) and C, 6 appear in f|T5l) . If M and T satisfy f|T6l) 
then C"C(M)(1 + T'^)T'^ < i, and hence the mapping is a strict contraction on 
Y. By Banach's fixed point theorem, ip has a unique fixed point n G Y, 

u{t) = U{t,0)(f)-i [ dsU{t,s)g{s,u{s)), (17) 



^We recall the results of Strichartz theorem, see [SI [12] and also [13]. For every G 
and every admissible pair (g,r), the function t C/(t,O)0 G L''(M^)) n C {R, L'^ {R'^ )) , 

such that 

l|(7(-,o)0|U,(K.L.) < CUh2, ycb e i2(KiV)^ 

where C is a constant that depends on q. Consider / C K such that e /. Let J C / such that 
G J, where " denotes the closure. Let (7,p) be an admissible pair, and f E L'' {I,Lf (R^)). 
Then, for all admissible pairs {q,r), the function 



t $/(t) = [ dsU{t, 
Jo 



s)f{s) e L«(/, L'-(M^)) n C( J, L2(M^)), 



such that 

\\^f\\Ll{I,L'-} < C||/||^^'(jj;,p'), 

where C is a constant independent of / and depends on q and 7 only. 
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for almost all t G {—T,T). Furthermore, since (p{u) G T, T], /f-^(M^)), u G 
C([-T, T],H\R'')). 

Step 2. Uniqueness. We now use the Strichartz estimates to prove uniqueness of 
the solution on the interval (— T, T). Suppose there exists u and v satisfying (IHl) 
on the interval (— T, T). Then 

w{t) := u{t) - v{t) = -i / ds U{t,s)[g{s,u{s)) - g{s,v{s))]. 

Jo 

Let ^ 

Wj{t) := -i j ds U{t,s)[gj{s,u{s)) - gj{s,v{s))], j = l,---,k. 
Jo 

It follows from Strichartz theorem that 

ft 

ds U{t,s)[gj{s,u{s)) - gj{s,v{s))]\\Lij{{-T,T).L''i) 



< C'\\gj{s,u{s)) - gj{s,v{s))\\^,>^^__^^^^^^,'^. 
Together with Assumption (B2), this implies 

II j ds U{t,s)[gj{s,u{s)) - gj{s,v{s))]\\Ly^^.T,T),L^^) 

< C'C{M){1 + T^)T '^"'^ \\u - t;|U'.((-T,T),L^.), 



Choose T < 1 small enough such that C'C{M)T '^'^ < i, j = 1, ■ ■ ■ , A;. Then 

\\'^j\\L''m~T,T),L''3(R")) < \\'Wj\\L''j {{-T,T),L''H^"))^ 

which implies Wj = on (— T, T), j = 1, ■ ■ ■ ,k. 

Step 3. Blow-up alternative. We now prove the blow-up alternative by contra- 
diction. We define 

T* := { sup T : a solution for ([8]) exists on [0,T*)}. (18) 

TeR+ 

Suppose that T* < oo and that 3M < oo and a sequence {tjjjeN C [0,T*) such 
that ti T* and ||n(ti)|| < M for all i G N. Choose j G N such that tj+T{M) > 
T*, where T[M) is the time scale over which Steps 1 and 2 holds. Applying the 
above analysis starting with u{tj) implies the existence of the solution of ([8]) to 
times tj + T{M), which contradicts f|T8|) . Therefore, ||m||hi — > oo if t y T*. We 
also define 

T* := { sup T: a solution for ([HD exists on (— T^,0]}. 
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Using a similar argument, one can show that ^ cxd as t \ — T*. Therefore, 

([8]) has a blow-up alternative. 

Note that it follows from ( flTl) . ( fT4l) and Strichartz theorem that 

«eLL((-r„T*),i^i'^(M^)), 

for all admissible pairs {a,b). 

Step 4- (ii) Charge conservation. To prove charge conservation, we use Assump- 
tion (B3) and the fact that u G H^. Using ([8]), we have 

II'^IIl^ ~ ^^""5 "^^t"") — (^""i —Am) + {iu, g{t, u)) = 0, 

and hence 

Mt)\\L- = U\\L'^, tG(-T„T*). (19) 

Step 5. (in) Continuous dependence. Consider the sequence {(f)n}nm such that 
0n — ^ in as n — s> 00, and let m„ be the maximal solution of (IHl) corresponding 
to the initial condition We claim that there exists T > 0, which depends on 
11011^1 only, such that m„, is well defined on [— T, T] for n large enough, and m„ — >• m 
in C([-T, T], LP(M^)) as n ^ 00, p G [0, ^), (p G [0, 00) = 1). Claim (Hi) 
follows by repeating this property to cover any compact subset of (— T=k,T*). 

Since 0„ — > in as n ^ 00, it follows that there exists no G N such that 
ll^nlli^i < 2||0||j:^i, Wn > riQ. It follows from Steps 1 and 2 above that there exists 
T = T(||0||j:/i) such that n„ and u are defined on [— T, T], for n > no, such that 

\\u\\l'^{{-T,T),H^) + sup \\Un\\L°°({-T,T),H^) < C''||0||/fi, 
n>no 

for some positive constant C. Note that charge conservation (ii) and 0„ — >■ in 
imply that 

Un^u in C{[-T,T],L\R^)). (20) 
Furthermore, it follows from f|T3l) that 

n(t)-n„(t) = f/(t,O)(0-0„) + (^(n)(t)-(^(u„)(t), tG [-T,T]. 

Using Strichartz estimates, 

||m — Mn.||L°°((-T,T)L2) + 11^ ~ Un\\Li({-T,T),L^) 

< C"||0 - 0„||^^i + C"C(M)(1 + T^)T'^x 

X {\\u - Un\\L°°{{-T,T),L^) + 11^ ~ Un\\ Li({-T,T),L^)) , 

where 6 = min,g|i ... k} -^-4^ > 0. Choosing T small enough (yet still depending 
on ||0||^^i) such that C'C{M){1 + T^)T^ < i, 

\\U - Un\\L^({-T,T),L^) + \\u - Un\\Li{{-T,T),L-) < 2C"||0 - 0n||f/i- (21) 
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It follows from (!20!) . (!2T]) and the Gagliardo-Nirenberg inequality that 

Ur^^u in C([-T,r],L^^(M^)), 

for all 'p G [2, -^l^^) (P ^ [0, oo) if = 1). This completes the proof of the 
proposition. 

□ 

We define the energy functional 

E{t,u) ■.= ^ J \Vu\^dx + G{t,u), (22) 

for u e H\R^). Note that since the nonlinearity and the potential depend on 
time, the energy is not conserved. We have the following proposition. 

Proposition 2. Suppose that Assumptions (B1)-(B4) are satisfied, and let u 
denote the solution of ^ given by Proposition Ui Then 

\E{t,um<\E{0,<P)\+TCm\L^), (23) 

for allt E [— T, T], where[—T,T] is a compact subset of {—T^,T*), and C{\\(f)\\L2) 
appears in Assumption (B4)- 

Proof. Since Assumptions (B1)-(B3) are satisfied, the results of Proposition [1] 
hold. We choose a finite T > such that T < min(T=i,, T*). We know from 
Proposition [1] that 

«GLL((-T„T*),iyi'^(M^)), 
for all admissible pairs {a,b). In particular, 

nGL''^((-T,T),iy^''-^(M^)), j = k, 

where the admissible pairs {qj,rj) appear in Assumption (B2). We note that by 
Mihlin's multiplier theorem, 1^*'^ = H'^''^ for 1 < p < oo and s an integer, see for 
example [221 . 

Since V commutes with the unitary propagator U corresponding to the free 
time evolution, and since the norm is invariant under unitary transformations, 
we have, using the Duhamel expansion of u given in (fT7|) . 

||VMt)||i. = ||Vf/(0,t)«(t)||i. 

= ||V0-z / ds U{0,s)Vg{s,u{sml, 
Jo 

t ft 

2 oT„,/v7j, / J„ TT^n „\V7„^„ , II / J„ rr/n „\V7„^„ „,^„NM|2 



||V0||i. -2Im(V0, / dsUiO,s)Vgis,uis))) + \\ / ds UiO, s)Vgis,uis))\\i, 
Jo Jo 

\\Vm, + 2lm[ ds {Vg{s,u{s)),U{s,0)V4>) + \\ [ ds U{0, s)Vg{s,u{sml, 
Jo Jo 
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Notice that 

|2 



rfs f/(0,s)V^(s,u(s))||i2 = 2Re / ds{V g{s,u{s)), / ds' U{s, s')V g{s\u{s'))) 

Jo Jo 

= 2Im / ds {'Vg{s,u{s)), —i / ds' U{s, s')^/ g{s' ,u{s'))), 
Jo Jo 

and hence 

II Vwlli^ = II V0||i. + 2 V Im r (V^?,(s, Vn(s)) 

= l|V</'||i2-2Vlm / (^?,(s,n(s)),An(s)), 

where the scalar product is well-defined using Assumption (B2) and duality on 
(Li((-T, T), ((-T, T), i/^'"^)) x (L°°((-T, T), H^)nL'i^ {{-T, T), iyi'^'O)), 

j = 1, ■ ■ ■ , A;, see for example [29]. Now, 

Im(^(t, u{t)), Au{t)) = limlm((l - eAy^g{t, u{t)), (1 - eA)~^An(t)) 

= limlm((l - eAy'git, u{t)), (1 - eAy\-zdtu{t) + g{t, ^(t)))) 
= limlm((l - eA)-^g{t,u{t)), - eA)-^9t?i(t)) 
= Re{g{t,u{t)),dtu{t)) 

= ^GitMt)) - {dtG){t,u{t)) 
dt 

for almost all t G (—T,T), where G appears in Assumption (B4). Therefore, 
l|Vn(t)||i2 = ||V</)||i2 -2G(t,u(t)) + 2G'(O,0) + 2 ! ds {dsG){s,u{s)). (24) 



Together with (1221) . Assumption (B4), and conservation of charge f[T^ . this im- 
plies 

|E(t,w(t))|<|E(O,0)|+TC'(||0|UO, 

for t G [— r, T] . The claim of the proposition follows by iterated application of 
this result to cover every compact subset [— T, T] C (— T*,T*). 

□ 

This is the main result of this section. 

Theorem 3. Suppose (B1)-(B5) hold. Then the solution u of ^ with initial 
condition cf) G H^{EJ^) is global in H^, i. e.,T* = T^, = oo, where T*, appear in 
Proposition [H 
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Proof. Since Propositions [T] and [2] follow from Assumptions (B1)-(B4), we only 
need to show that t G [0,T*) is finite if T* < oo, which, together with 

the blow-up alternative, implies a contradiction. Suppose T* < oo. Assumption 
(B5), (Hn]) and ((231) imply that 

^h(^)llHi = ^(h(^)lli^ + l|Vn(.)||i.) 

<^\\u{s)\\l. + \E{s,u{s))\ + \G{s,u{s))\ 

< Imh + \m 0)1 + T*cmL^) + ^-^ms)\\i. + cdi^iuo, 

for all s G [0, T*), and hence 

< oo, 

for finite T*, which contradicts the blow-up alternative. The case of T* is proven 
similarly. 

□ 

Remark 3. One can directly verify that Assumptions (Al) and (A2) in Suhsec- 
tion \2.S\ imply that Assumptions (B1)-(B5) are satisfied. 

Remark 4. Assumptions (B1)-(B5) are satisfied if 

git, = Vit, >(■) + /(t, -, u{-)) + iWit) ^ |n|2)(>(.), u e H\ (25) 

where V, f and W satisfy the following; see f27^ . 

The potential V is real valued such that V,W G C(M,Lp(R^) + L~(M^)), 
with p > f,p > 1, and dtV G L°°(M, L°°(M^)). Let giit,u){-) = V{t,-)u{-). 
Assumption (B2) follows by Sobolev's embedding theorem with ri = ^zi^P = 1) 
and Assumption (B3) follows trivially since V is real valued. 

The local nonlinearity / : M x x [0, oo) — > M such that f{t,x,u) is con- 
tinuous in t, measurable in x and continuous in u, and 0) = OVt G M and 
almost every x G M^. If N > 2,3 a positive constant C and a G [0, ^^y^^) such 
that 

\f{t,x,u) - f{t,x,v)\ < C(l + \ur + \vn\u-vi 

uniformly in t G M., for almost all If N = 1, then for every 

M > 0, 3L{M) such that 

\f{t,x,u)- f{t,x,v)\<L{M)\u-v\, 

uniformly m t G M, for almost all x G M.'^,u,v G M, |u| + |f| < M. The local 
nonlinearity f is extended to M. x x C by defining 

f{t,x,z) :=±f{t,x,\z\), V2GC\{0}, 
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t G M and for almost all x G M . We define 

g2it,u){-):= f{t,;u), teR,uEH\ (26) 

and 

r rH 

G2{t,u) = J dx J drf{t,x,r), ueH^, t G M, 

Assumption (B2) follows from Sobolev's embedding theorem ^ L""''^(]R^), r2 = 
a + 2. Assumption (B3) follows from [2^) . 

For the nonlocal (Hartree type) nonlinearity {W -k \u\'^)u, W is a real valued 
function, W -.RxR^ ^R, such that W G L'?(M^), g > f , g > 1. We make the 
identification g^{t,u) = {W{t)'k\u\^)u, t G M, m G H^, andG'j,{t,u) = \j dxiW-k 
\u\^)\u\^. Assumption (B2) follows from Young's and Sobolev's inequalities, with 
= 2^3^. Note that the analysis can be directly extended to the case when f is a 
finite sum of local and nonlocal nonlinearities. 

Assumptions (B4) and (B5) follow by the application of Holder's and the 
Gagliardo-Nirenberg inequality if the above holds together with the following. There 
exists 5 G [0, jf) and C > constant such that 

drf{t,x,r) < G\uWl + \u\^), 
uniformly in t eR, and 

dt /' 'dr/(^,x,r)<A(t)|n|^ 

for t eR, almost all x e R^ , and u e H^, where A{t) G L°°(M). Moreover, W 
in the nonlocal nonlinearity is spherically symmetric such that 

:= max(0, W) G L°^(M, L'^(M^) + L°°(M^)), 

dtW+ G L°°(M,L°°(M^)). 

Remark 5. There are several ways in which one may relax condition (B4)- The 
above result also holds for l{t) appearing in (B4) in L^{R), 1 < p < oo. If l{t) 
appearing in assumption (B4) is in L^{R), then one can easily show that the 
norm of the solution of the nonautonomous nonlinear Schrddinger equation ^ 
is finite for all times t E R. In this case, one may extend the above analysis to 
potentials which grow at infinity, say quadratically; see IWf for a discussion in 
the time-independent case. 

On the other hand, if l{t) G L^^^{R), 1 < p < oo, the upper bound on the 
energy functional, / f^) in Proposition is replaced with 

m^um < |E(o,0)| + rc"(r)C'(||0|U.), 
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where C'{T) > is finite for finite T. 

Furthermore, if one replaces Assumption (B4) with 

\dtG{t,u)\ < \\u\\Lr>l{t), 




we have 



\E{t,um<\E{0,<P)\ + C{T), 



where C{T) > is finite for finite T. Global well-posedness follows like before. 

Remark 6. There are relatively few rigorous results on nonautonomous nonlin- 
ear Schodinger equations. The Cauchy problem for a nonautonomous nonlinear 
Schrddinger equation that is obtained by applying a pseudo-conformal transfor- 
mation to a nonlinear Schrddinger equation with local nonlinearity is studied in 
l31^ . Moreover, well-posedness of a nonlinear Schrddinger equation with a lo- 
cal nonlinearity whose coefficient is time-periodic was studied in l3^ . where the 
instability of the ground state due to periodic modulation of the nonlinearity is in- 
vestigated. Furthermore, (endpoint) Strichartz estimates were obtained in JS^ for 
time- dependent potentials which are small and concentrated in frequency space. 
Time- dependent potentials also arise in the analysis of charge transfer models , 
f34^ , and scattering of multisolitons, fS^ . where Strichartz estimates were applied 
in order to study the asymptotic stability of multisolitons. We also mention in 
the linear case the dispersive estimates in 13^ for time dependent potentials that 
decay in time, and the analysis in [37^ on the slow growth of Sobolev norms for 
the linear Schrddinger equation with (quasi-) periodic potentials. 

4 Properties of the nonlinear Schodinger equa- 
tion 

In this section, we recall some properties of the nonlinear Schrodinger equation 
([1]) and the soliton manifold, see for example [27]. We will use these properties 
in the following sections. 

4.1 Symplectic, Hamiltonian and Variational structure 

The space H^{M.^,C) = if^(]R^,M^) as a real space, and it has a real inner 
product (Riemannian metric) 




(27) 
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for u, V It is equipped with a symplectic form 

uj{u,v):=lm / dxuv={u,iv). (28) 



The Hamihonian functional corresponding to the nonhnear Schrodinger equation 
© is 

Hvii^) ■■=\j (I V7/;|2 + V\^'')dx - Fii;). (29) 

Using the correspondence 

H\R^,C) < — > H\R^ ,R) ® H\R^ ,R) 
ip < — > {Keip, Imip) 

i — > J, 

where J := (^^i complex structure on if^(M^,]R^), the nonlinear 

Schrodinger equation can be written as 

In the following, we denote H\R^) as either H\R^ , C) or H\R^ , M^). We note 
that since the external potential is time-dependent, the Hamiltonian functional 
Hv defined in ( l29l) is nonautonomous, and there is no conservation of energy. 
Still, Hy is invariant under global gauge transformations, 

and the associated conserved Noether charge is the "mass" 

iV(^) ■.= ljdx (30) 

Orbital stability (Assumption (A5), Subsection 12. ip implies that appearing 
in Assumption (A4) is a local minimizer of Hv=o{'ip) restricted to the balls Bm '■= 
{ip G : N{iIj) = m}, for m > 0; see [T71 [T8] . They are critical points of the 
functional 

■■=lldx (I V^p + fi\^\') - F(^), (31) 

where /i = fi{m) is a Lagrange multipher. 
3 The tangent space TH^ = H^. 
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4.2 Soliton Manifold 

When V = 0, Assumption (A3) implies that the nonhnear Schrodinger equa- 
tion ([T]) is invariant under spatial translations, time translations, gauge trans- 
formations, spatial rotations and Galilean transformations. The corresponding 
conserved quantities are the field momentum, energy, mass, angular momentum, 
and center of mass motion, 

^j\V^\^-F{ij), J^{xA-tV)^, j i,{x + 2itV)ip. 

When 7^ 0, the above quantities are generally no more conserved. In 
particular, the rate of change of energy is 

dtHvW = ^j dxdtV\^\^ (32) 

and the rate of change of momentum is 

dt{ij,-zVij) = -{i>,{VV)ij). (33) 

Formally, Eq. (!32l) follows from ([1]) and (129|) . while (!33|) . which is a statement of 
Ehrenfest's Theorem, follows from ([1]). We refer the reader to Appendix [7] for a 
proof of (132]) and (|33|). 

We introduce the combined transformation Tavy- 

i^av, ■■= Ta^y^P = e*(^''-("-'^)+^V(a^ - a), (34) 

where v,a E and 7 G [0, 27r). We define the soliton manifold as 

Ms := {Z,a.yV,. : a, i;, 7, /i G X X [0, 27r) x /}, (35) 

where / appears in Assumption (A4). If /'(O) = 0, where / appears in (Q, then 
/ C M"^. The tangent space to the soliton manifold Ais at 77^ G Aig is given by 

Tri^Ms = span{et, Cg, e^, e J, (36) 

where 



eb 



Remark 7. In the case of pure local nonlinearities, f{ip) = Xlipl'^ip, and V = 0, 
the nonlinear Schrodinger equation (QP is invariant under scaling, 

■ '^{x, t) fJ'~^i^{VJ^x, fit). 

In this case, one can define the generalized transformation Tavy fi '■= TavyT^. Fur- 
thermore, Cf, = df,T^r]\f,=i = |(f + x ■ V)?7i. 
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Remark 8. Note that the solitary wave solution rj^ breaks the translation and 
gauge symmetries of the nonlinear Schrddinger equation, which leads to associ- 
ated zero modes of the Hessian defined in Assumption (A6), Suhsection \2.1\ 
Differentiating £'^{T^T^ri^) = with respect to a and 7 and setting the latter two 
to zero gives 

C^Ct = , C^Cg = , 

while a direct computation gives 

Remark 9. The soliton manifold M.^ inherits a symplectic structure from {H^, u). 
For a = {a, v, 7, /i} G x x [0, 2%) x J, 



where P^ is the orthogonal projection onto Tj^^M-g. It can he shown that VL^^ 
is invertible if d^m^jj) > 0, where the mass m(/i) = ^ f dx rf^, see Explicitly, 



( ° 


— m(/i)l 





\ 


m(/i)l 




















m'{fj,) 


V 





—m'{fi) 


/ 



2N+2 



(37) 

where ek,k = 1, ■ ■ ■ , 2N + 2, are basis vectors of %i^M.s (^nd 1 is the N x N 
identity matrix. Note that is related to by a similarity transformation. 

Remark 10. It has been noticed in ^ that in one dimension, the action of 
the combined transformation on elements of the soliton manifold has a group 
structure, the Heisenberg group . In N -dimensions, the group corresponds to 
l_|2Af+i^ t/ie Heissenberg group in 2N + 1, which is given by 

(a,w,7) ■ (a ,w',7') = (a",w",7"), 

where a" = a + a', v" = v + v' , ;7" = 7' + 7 + \va' . Note that the Heisenberg 
group is a central extension of the additive group. 



4.3 Skew-Orthogonal Decomposition 

Consider the manifold A^', = {r/^, a G Sq}, Sq = x x [0, 27r) x Jq, where 
/o C I\dl is bounded. We define the 5 neighbourhood of M'^ in as 

Us:={^ljeH\ inf llV^-r/.ll <5}- 

Then, for 5 small enough and for all ip G t/^, there exists a unique (t{iI)) G 
C^iUs.T.) such that 

^^(^ - ^/aW, e) = ?7aW, J~^e) = 0, 

for all e G %i^^^^Ais- For a proof of this statement, we refer the reader to [1], see 
also [391. 
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5 Proof of the main result 



The proof of the main result is based on an extension of the analysis in |4] , except 
that one needs to keep track of additional terms due to the time-dependence of 
the potential. Formally, the proof boils down to decomposing the solution into 
a component which belongs to the soliton manifold plus a fluctuation which 
is skew-orthogonal to the soliton manifold. The dynamics of the component 
belonging to the soliton manifold is effectively determined by the restriction of 
the Hamiltonian flow generated by the nonlinear Schrodinger equation to the 
soliton manifold, while the norm of the fluctuation is controlled using an 
approximate Lyaponuv functional. An additional feature of our analysis is an 
iteration scheme which gives a much longer time scale, 0{\ logh\/h), over which 
one can control the fluctuation, compared to 0{h~^) in |4]. 



5.1 Reparametrized equations of motion 

Suppose ip satisfies the initial value problem ([1]) such that ip E Us C , where Us 
appears in Subsection 14.31 By the skew-orthogonal decomposition, there exists a 
unique a = ct(^) = {a, w, 7, e S = x x [0, 2tt) x / and w' G such 
that 

ip = r]^ + w' , 
w' ± J-%^Ms. 

Let 

u ■■= T'l^il) = r]i, + w, (38) 
where w = T'J^r]^. We introduce the anti-selfadjoint operators 

Lj = dxj, Ln^j = —Jxj, j = 1, ■ ■ ■ , N, 

L2N+1 = —J, L2N+2 = d^, (39) 

with coefficients 

a J = dtaj - Vj, aN+j = ~dtVj - d^V^t, a), j = !,■■■ ,N, 

«27v+i = /i - ^v'^ + ^dta ■ V - V{t, a) - dt'j, a2N+2 = -dtfJ'- (40) 
We denote by 

|a| := sup \ai\, (41) 

je{l,--- ,2N+2} 

and C{a,w,h) := + h"^ + We have the following proposition. 
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Proposition 3. Suppose that ip G satisfy ^ such that ip{t) G Us for t G 
[0,T]. Then the parameter a = {a,v,-f,fi} and w G H^, as given above, satisfy 

dtaj = Vj + 0{C{a,w, h)) 

dtvj = -2d,y{t, a) + 0(C(a, w, h)) 

dt'j = fi- ii;^ + ^dta ■ V - V{t, a) + 0{C{a, w, h)) 

dtfi = 0{C{a, w, h)), 
/ortG(0,T),j = l,---iV. 
Proof. We proceed in three steps. 

Step 1. Equations of motion in center of mass reference frame. 

We first find the equation of motion for u = T^l^ijj = e~2^^''^^"'^ip(x + a). 
Differentiating u with respect to t and using ([T]), and the fact that 

g-f («.x+7)^^^^ + a) = An + ■ Vu - — 

e-f(--+^)/(V;(x + a)) = /(«), 

we have 

2Ar+l 

= - J((-A + - f{u)) + ajLjU + JTZyu (42) 

where 

TZv := V{t, x + a)- V{t, a) - VV{t, a) ■ x, 
and Lj and aj are as defined in and (HQ]) respectively. In other words, 

2N+1 

dtu = J£'^{u) + ajLjU + JTZyu, 
i=i 

where £^ appears in fl3Tl) . 

Step 2. Reparametrized equations of motion. 

We now use P2|) and the skew-orthogonal decomposition to find equations for 
the parameters a = {a,v,'-f, fi} and w. Recall that 

^',iv,) = 0, 

which implies 

where C, = (-A+^-f (r/^)) = S;{r^,) and N,{w) = f{^^ + w)-f{r^,)-nr^,)w. 
Substituting back in fH2l) implies 

2N+1 2N+2 

dtw = {JCf, + ^ ajLj + JTZv)w + Nfj,{w) + J ^ (^jLjf]^ + JTZyr]^. 
i=i i=i 
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We know that {Jz,w) = for all z G T^Ais- It follows that 

dt{Jz, w) = dtfi {Jdf,z, w) + {Jz, dtw) = 0. 

Therefore, 

dtfi {Jdf,z,w) = -{Jz,dtw) 

= -{Jz, JC^w) - {Jz, JTZviVf^ + w)) - {Jz, Nf,{w))- 

2N+1 2N+2 

-{Jz,Y^ ajLjw) - {Jz,Y^ ajLj7]f,). 
i=i i=i 

It follows from Remark ([H]), Section HI that {Jz, JC^w) = 0. Together with 
[L,,J] = 0, L* = -L„ j = l,--- ,2N + 2, 

this implies 

2N+2 2N+2 

^ aj{JLjZ,w) = -{z,nviw + 7]^) + Nf,{w)) + {Jz,ajLjr]^). (43) 
i=i i=i 

Choosing z = e^, where e^, /c G {1, ■ ■ ■ , 2N + 2}, is a basis vector of T^jAis gives 

2N+2 2N+2 

= (cfc, iV^l^) + T^viw + ?7/,)) + ^ aj{Ljek, Jw), 

where appears in ( !37j) . Replacing the definition of aj, appearing in (l40l) . in 
l3l) . and using the fact that 



{xjrif^jnvr]^) = 0, rii,{x) = Vti{\x\), {Jr]f„TZvVti) = 0, 

gives 

^ 27V+2 

<9tafc = Wfc H ^{{xkVfj,, JN^{w) + JTZyw) + ^ aj{LjXk'niM,u))), (44) 

(9ti;fc = -2(9a;^\^(t,a)H i-^{{dx^ri^,N^{w)+TZvw)- ^ aj{Ljdx^r]^, Jw) + {dx^r]^, 

(45) 

9*7 =/i - ^t;^ + ^SfO ■ - V(t, a)- 

^ 2Af+2 

T7—:{.{^^lV^Ji^N^,{w) +'Rv{w + ri^)) - ^ aj{Ljdf,ri^, Jw)), (46) 
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2N+2 



(?7^, JNf,{w) + JTZyw) - ^ aj{Ljr]f„w). 



(47) 



The claim directly follows from Assumptions (Al) and (A2), which imply that 



5.2 Control of the fluctuation 

In this subsection, we use an approximate Lyaponuv functional to obtain an 
explicit control on and This approach dates back to [HI [20], and has 

been used in [H |5] . We define the Lyapunov functional 



where is defined in ( 13T|) . We proceed by estimating upper and lower bounds 
for C^{u, r]^), where u = t]^ + w appears in (1551) . 

5.2.1 An upper bound for the Lyapunov functional 

Lemma 1. Suppose ip satisfies (QP such that ip{t) G Us, t G [0,T], for some 
(5 > 0, and let u, w, Tj^ he as defined in Subsection 15. ii Then there exists a 
constant c independent of h such that 



WRvekWL^ = 0{h^), \\N,{w)\\h^ < C\\w\\hi, 
for < 1; see Remark[T], Subsection 12. 1[ 



□ 



C^{u, v) := S^{u) - S^{v), u,ve H\R^), 



(48) 



C^,{u,7]^,) < ct{h'^\\w\\m 



+ 



{\a\ + h)\\w\\lr), 



(49) 



where \a\ appears in ( [7?! ). 
Proof. Recall that 




By translational symmetry, 




where V-a{,x) = V{x — a). Furthermore, 



Hv{T-\^) = HvW + l{\v' + /i) ll^lli. 



v ■ {iip, Vip) + 
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and hence 

S,{u) = Hviij) + ^i^v' + ^) llV'lli. - ■ {tij, dxV\^\ (50) 

We have the following relationships regarding the rate of change of field energy 
and momenta. 

dt^j dxV\^lj\^ = {VVii^.Vilj) + ^{ilj,dtVij), (51) 
dtHv{'ilj) = ^{ilj.dtV^) (52) 

and Ehrenfest's theorem 

ai(#,V^) = -(^,V\/V')- (53) 

Using ([T]), the above three statements flM]) . flS^ and flS51) are formally verified. 
Rigorously, one can verify them using a regularization scheme and a limiting pro- 
cedure, and we refer the reader to the Appendix for a proof of these statements; 
see also the proof of Proposition [21 Section [31 Furthermore, it follows from gauge 
invariance of ([T]) that the charge is conserved, 

dt\mL^ = Q, (54) 

see Theorem [3l in Section [3l Differentiating (l50l) with respect to t and using 
(Eiai) gives 

X 1 1 

-{VVt^,Vij)-~{ij,dtVij) 

= ^(^^ + dtf^)Mh - \dtv ■ VV^) + \v ■ (^, VVij) - (VVz^, VV^) 
= \dtii\\u\\l, - {]^{i{dtV + 2WVa)u,Wu), (55) 

where we have used u{x) = e**-2^'^^'*''*'?/'(x + a) and translation invariance of the 
integral in the last line. Furthermore, it follows from (!3T|) that 

which, together with (HSjl and ([S5D, implies 

dtC^{u,7],) = IdtfiiMl. - WvXl^) - {t{^dtv + VVa)u,Vu). (56) 
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We now estimate both terms in (156!) . Since {iz,w) = for all z G T^Ais, it 
follows from the skew-orthogonal decomposition (Subsection 14.31) that 

Mlh - WvJh = \\w\\l2, 

and hence 

^5.Mlklli.-hMlliO = 0(|a|||^||i.). (57) 

To estimate the second term, we replace u = r]^ + w, and use the fact that 
(icg^w) = (iV?7^,w) = 0, and {ihri^,Vri^) = for all real h G L°°(]R^). We have 

{ii^dtv + VK)m, Vm) = {iVVaW, Vr/^) + {^VVar]^., Vw;) + {i{^dtv + VVa)w, Ww). 
Adding and subtracting W(t, a) ■ {iw, Vr/^) = W{t, a) ■ {irj^, Vw) = gives 

{^{l^tv + VVa)u, Vn) ={^dtv + VV{a)){iw, Vw) + ((VK - VV{a))zw, Vw) 

+ ((VK - VVia))tr]^, Vw) + ((VK - VVia))tw, Vr/^). 

The first term of the above equation is of order 0(|q!| ||w||^i), while Assumption 
(A2) implies that the second term is of order 0(/i||w||^i). Assumptions (A2) and 
(A4) imply that the third and forth terms are of order 0(/;,^||u;||hi). Hence the 
claim of the lemma. □ 

5.2.2 A lower bound for the Lyapunov functional 

In this subsection, we estimate a lower bound for C^{u, rj^). Let 

:= {w G H\R'') : {w, rh) = 0, Vz G %,^Ms}. 
It follows from the coercivity property of that there exists a positive constant 

p := inf {w,CM > 0- (58) 

We refer the reader to the Appendix D in for a proof of this statement. We 
have the following result. 

Lemma 2. Suppose tp satisfies (QP such that tpit) G Us, t G [0,T], for some 
6 > 0, and let u,w,ri^ he as defined in Suhsection \5.1[ Then there exists positive 
constants p and c independent of h such that, for \\w\\h'^ < 1, 

\C^{u,r]f,)\ > ^\\w\\hi - c||u;||f^i, (59) 

where p appears in ^5^. 
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Proof. We first expand S^{u) around 77^, which is a critical point of Sf^. 

+ w) = £M + + Rf\w), (60) 

where 

R^;^\w) = F{r]^ + w)- F{7],) - {F'{7]^),w) - ^{F"{r]^)w,w). 

It follows from Assumption (Al) that 

\R^^\w)\ < c\\w\\l,, 

for llu'll/fi < 1, where c > is independent of t G M; see Remark [T] in Subsection 
12.11 Furthermore, the coercivity property (158|) implies 

{w,C^w) > p\\w\\hi, 

and hence 

K{u,ri^)\ = \Sf,{u) - Sf,{rif,)\ > - c||w||^i, 

for < 1- 

□ 

5.2.3 Upper bound on the fluctuation 

In this subsection, we use the upper and lower bounds on the Lyapunov functional 
to obtain an upper bound on 

Proposition 4. Suppose (Al)-(A7) hold, and let ip satisfy ([2]j, and u,7]p^,w as 
above. For h <^ 1, choose T G M"*" such that ip{t) G Us, t G [0,T], where 
Us appears in Subsection \4.3\ Fix e G (0,1), and choose to ^ [0,7"] such that 
||w(to)||^i < Then, fork small enough, there exist absolute constants Ci> 1 
and C2 > 0, which are independent of h,e and t, such that 

sup |k(t)||^i <Ci(/l2+||w;(to)|||0) 

t£[to,to+T] 

sup \a^{t)\<Ci{h^+\\wito)\\jji), 

ie{l,--- ,2Af+2} 

where r = C2/h. 

Note that the conditions of Proposition |4] are satisfied for to = 0. 
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Proof. It follows from Lemma [T] that, for t >to, 

\C^{Vf^+wit),T]^)\ < |C^(?7^+tf;(to),?7^)|+c(t-to)(|a||k(t)||^i+/ilu;(t)||//i+/i||tf 

Furthermore, expanding S^{r]n + w{to)) around 7]^ and using Assumption (Al) 
gives the upper bound 

\C^,iv^^ + '^^ito),r]^)\ < c'\\w{tQ)\\jji, for \\w{to)\\m < 1, 

where c' is a constant independent of h,e, see Remark [H Subsection 12. 1[ There- 
fore, 

\C^{r]^ + w{t),r]^)\ < C\\w{to)\\li + C{t - to){h^w{t)\\H^ + (|a| + h)\\w{t)\\]j^), 

for some constant C independent of h and e. Together with Lemma [2], it follows 
that 
1 



where p appears in (!58|) . Equivalently, there exists a positive constant C inde- 
pendent of h and e such that 

C\\w{t)\\H^ < \\w{t,)\\l, + {t- to){h''\\w{t)\\H^ + (|a| + h)\\w{t)\\l,) + \\w{t)\\l,. 
^ - ^0 < 2(hq^ =: 

C\\w{t)\\l. < \\w{t,)\\l. + ^h\\wmm + f IkWII?,. + Wwm'm- 
Using the fact that 



h\\w{t)\U. <\h' + \\\w{t)\\l,, 



we have 

+ jh^ - j\\w{t)\\l, + > 0. 

Let I/O := \\w{to)\\Hi, y := supt^^t^^^^^^^ \\w{t)\\Hi, and f{y) = - ^y^ + yi + jh^. 
For /i <C 1, the function intersects the x-axis in a point such that < ci{h'^ + 
yg), where ci is a positive constant independent of h and e; see Figure 1. It follows 
that for yo < y^:, y < y^,, for t G [to, to + r]. Substituting back in (14411471) and using 
(HOj) and dS]) gives 

|a| < C2(/i^ + i/o), 

for some positive constant C2 which is independent of h and e. It follows that for 
h small enough, there exists positive constants Ci and C2 which are independent 
of h and e, such that 

sup \\w{t)\\l-, < Ci{h^ + \\w{to)\\l,) 
sup \aiit)\<Ci{h^+\\w{to)\\l,). 

te[to^to + -T-l 
je{i,--- ,2Af+2} 
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This concludes the proof. 




5.3 Proof of Theorem H 

In this subsection, we prove Theorem [2]by iterating the apphcation of Proposition 
m and using the resuh of Proposition [3d 

Proof of Theorem\^ Fix e G (0, 1). For h small enough, let T* be the maximal 
time for which the skew-orthogonal decomposition is possible; see Subsection 14.31 
Consider the interval 

[0, T] = [to, tl] U [ti, ts] U ■ ■ ■ U [tn-U tn] C [0, T*], 

such that 

C2 

= to <ti < ■ ■ ■ <tn = T, {U+i -ti) < —, ? = 0, ■ ■ ■ , n - 1, 

n 

where C2 appears is Proposition HI We will choose n E N depending on h and e 
later. Let 

\a\i := sup \aj\, 

je{l,-,2N+2} 

yj := sup \\w{t)\\Hi, j = 0, ■■■,n-l. 

t&[tj,tj+i] 

Note that yo < h and |a|o < Ch, for some constant C independent of h and e. 
Iterating the application of Proposition H] we have 

n 

yl<(J2ci)h'<cr'h' 
i=i 

\a\n < CC^-^^h\ 



''This iteration scheme is similar in spirit to the one used in [38]; see also |7]. 



32 



We choose n such that C^'^^h? < This imphes 



n + 1 < -e 



\ogh 



Therefore, for t G [0, e 



log Ci h 



\w{t)\\j,^<h- 



2-e 



The claim of the theorem follows by the application of Proposition [3l 



□ 



6 Two physical applications 

In this section, we sketch two physical applications of our analysis. Through- 
out our discussion, we suppose Assumptions (Al)-(A7) are satisfied, so that the 
results of Theorem [2], Section \2.2\ hold. 

6.1 Adiabatic transportation of solitons 

We discuss adiabatic transportation of solitons in time-dependent confining po- 
tentials. Suppose the external potential is locally harmonic and decaying when 
||a;|| oo such that 



where s G M and ^ ^ 1. Suppose the soliton center of mass is initially at a: = 0. 



Wh{t,x) = h^ujl{x- St), \\x-st\\<e 




Making the change of variables 



a = a — st 



V = V — s 



the above equations become 



d^h = v + 0{h^), 

dtv = -ujlh'^a + 0(0), 
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whose solution is of the form 



( Q j — g 3 

aU) = a(0) cos(huJot) + sm(hujot) + Oih"^). 



For a(0) = u(0) = 0, 



S 3 

ait) = st — - — sinihuoot) + 0(h^^ 
hujn 



for t < C\\ogh\/h. If U = 0(/i°),a > |, a{t) = st + 0{hl) ,i.e. the soliton 
is transported adiabatically, up to error terms due to radiation damping and 
oscillations. 



6.2 Mathieu instability due to time-periodic perturbation 
of trapped solitons 

As a second application of our analysis, we discuss the onset of Mathieu instability 
of initially trapped solitons due to time periodic perturbations. We note that 
our analysis can be easily generalized to more general time-periodic potentials, 
[101 EH], and quasi-periodic potentials, [H]. 

Suppose the external potential is locally harmonic and decaying at spatial 
infinity such that 

Wh{t,x) = h^culil + 5cos{ujt))x, \\x\\ <9 

where 6' 1. We claim that for special values of ujq/uj the system exhibits 
Mathieu instability, in the sense that the minimum of the potential becomes 
unstable under small perturbations. Choosing e = | in Theorem [21 we have, for 
9 large enough, and time t < C\ logh\/h, 

3 

dta = v + 0{h^), 

dtv = -ujI{1 + 6cos{ujt))h'^a + 0{h^). 

Note that this is nothing but Mathieu's equation, plus error terms due to radiation 
damping. For more general periodic potentials, one obtains Hill's equation, |40j . 
The system exhibits parametric resonance if huo = n^, n = 1, 2, ■ • • , 5 > 0, see 
for example [39], Chapter 5. Although the nonlinear term acts as friction, the 
minimum of the potential is unstable for 6 large enough, with /i2 ^ 5 <C 1. 
We note that the analysis is restricted to the instability of the minimum of the 

3 

potential. When ||a|| grows to 0{h~i) due to the instability, neglecting the 
nonlinear terms in the equations of motion is no more justified, and one needs 
to analyse the full nonlinear problem, see for example [12]. Instead, depending 
on the form of the potential, one can perceive different scenareos: The center of 
mass of the soliton oscillates chaotically, or there is a bifurcation and new points 
of stability. Further discussion of this application will appear elsewhere. 
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7 Appendix 



Rate of change of field energy and momenta, Section [5] 

Proof of ( f57]) . Differentiating Vip) with respect to t and using ([T]), Assumption 
(Al), and the fact that tp e i7^(M^), we have 

= iij, dtViP) + {-Aij + Vij- fi^p), tV^P) 

+ {t^,v{-AiJ + v^-f{m 

where we have used integration by parts in the last step. 

Proof of [5^] . The proof of (152|) follows directly from (!24|) in the proof of Propo- 
sition [21 Section [3l with the identification 

g{t,u) = V{t)u- f{u). 

Proof of [53\] . We use a regularization scheme similar to the one used in Propo- 
sition El Section [31 Let /, := (1 - eA)^\ 

dtiiip., Vip) = dt \im{IJ'ip, hVip) 

= \im{{hidtiJ,hVij) + {hi^^hVdtij)] 

= lim{(/,(-AV' + Vi,- /(^)),/,V^) - {h^,hV{-A^ + V^- f{'4j)))} 

= {ij,VVij)-{ij,V{Vij)) 
= -{ij,VV,p). 
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